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Inclusive transverse momentum distributions of charged hadrons within 0.2<j>t<6.0 GeV/c have been mea- 
sured over a broad range of centrality for Au+Au collisions at % /snn=130 GeV. Hadron yields are suppressed 
at high pr in central collisions relative to peripheral collisions and to a nucleon-nucleon reference scaled for 
collision geometry. Peripheral collisions are not suppressed relative to the nucleon-nucleon reference. The sup- 
pression varies continuously at intermediate centralities. The results indicate significant nuclear medium effects 
on high pr hadron production in heavy ion collisions at high energy. 

PACS numbers: 25.75.-q 



QCD predicts a phase transition at high energy density 
from hadronic matter to a deconfined Quark-Gluon Plasma 
(QGP) [jjj]. This transition may be studied in the laboratory 
through the collision of heavy ions at ultrarelativistic ener- 
gies. Partons propagating in a medium lose energy through 
gluon bremsstrahlung [[| |J], with the magnitude of the en- 
ergy loss predicted to depend strongly on the gluon density of 
the medium. Measurement of partonic energy loss therefore 
provides a unique probe of the density of the medium. 

Analysis of deep inelastic scattering and Drell- Yan pair pro- 
duction using nuclear targets indicates that the energy loss 
in cold nuclear matter is 0.2 — 0.5 GeV/fm for quarks with 
energy greater than 10 GeV [j^, Hard scattering of par- 
tons in nuclear collisions occurs early in the evolution of the 
extended system, thereby probing the phase of highest den- 
sity. Energy loss softens the fragmentation of jets, leading 
to the suppression of high transverse momentum (high pr) 
hadrons in the final state |Q] . The PHENIX collaboration has 
reported the suppression of charged hadron and ir° produc- 
tion at high pt in central Au+Au collisions at center-of-mass 
energy per nucleon pair - v /snn=130 GeV, relative both to ref- 
erence data from nucleon-nucleon (NN) collisions and to pe- 
ripheral Au+Au collisions The suppression is in quali- 
tative agreement with predictions of partonic energy loss in 
dense matter, though quantitative conclusions require the un- 
derstanding of other nuclear effects [^J . 

This Letter presents a measurement of the inclusive charged 
hadron yield (h+ + h~)/2 within 0.2<p T <6.0 GeV/c, mea- 
sured for a broad range of centrality in Au+Au collisions at 



^/si^BO GeV by the STAR collaboration at RHIC. Sup- 
pression of charged hadron production at high pr in cen- 
tral collisions is observed, in qualitative agreement with the 
PHENIX measurement [g]. The high precision and wide kine- 
matic and centrality coverage of the data presented here per- 
mit a detailed study of nuclear medium effects on hadron pro- 
duction from the soft to the hard scattering regime. 

For comparison of spectra from nuclear collisions to an NN 
reference, the nuclear modification factor is defined as 



Raa{pt) 



d 2 N AA /dp T dri 
T AA -d 2 a NN /dp T di 1 



(1) 



where Taa= (Nbin)/^^ account s for the collision geometry, 
averaged over the event centrality class. (Nt>m)> the equiv- 
alent number of binary NN collisions, is calculated using a 
Glauber model. Raa(pt) is less than unity at low pt [§]. In 
contrast, the yield for hard processes scales as (Nbi n ) in the 
absence of nuclear medium effects (Raa(pt)=1), and effects 
of the medium may be measured at high px by the deviation 
of Raa{pt) from unity. In addition to final state energy loss, 
Raa{pt) may be affected by initial state multiple scattering 
|Q, [l(J and transverse flow, both of which will enhance hadron 
production at high px, and by shadowing of nuclear parton 
distributions. At significantly lower y/s than the present study, 
enhancement of hadron production at high pt has been ob- 
served in p-nucleus [ 1 1 1 and a — a Jl2| ] collisions, as well as 
central collisions of heavy nuclei Jl3|]. Current estimates of 
shadowing effects at RHIC energies contain large uncertain- 
ties \fi[ p~ 
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TABLE I: Typical multiplicative correction factors and systematic 
uncertainties, applied to the yields for peripheral and central colli- 
sions. 



centrality 


p T =2 GeV/c 


PT =5.5 GeV/c 


60-80% 


0-5% 


60-80% 


0-5% 


Tracking 
Background 
Pt resolution 


1.16±0.08 
0.95±0.03 
1 00 +0 


1.59±0.16 
0.90±0.05 
0.99 ±0.01 


1.14±0.23 
0.95±0.05 
0.92±0.05 


1.49±0.30 
0.88±0.12 
0.76±0.10 



The STAR detector is described in [p"5j]. Data collection 
utilized both a minimum bias trigger and a trigger selecting 
the 10% most central events. Charged particle tracks were 
detected in the Time Projection Chamber (TPC), with mo- 
mentum determined from their curvature in a 0.25 T mag- 
netic field. After event selection cuts, the central data set 
contained 320,000 events, while the minimum bias data set 
contained 240,000 events. The measured minimum bias dis- 
tribution, corrected for vertex finding efficiency at low multi- 
plicity, corresponds to 94±2% of the Au+Au geometric cross 

tAuAu 



section <7g e "m , assumed to be 7.2 barns [|1 

Centrality selection is based on the primary charged parti- 
cle multiplicity N c h within the pseudorapidity range |?y|<0.5. 
The most central bin is 0-5% of cr^,^, while the most pe- 
ripheral bin is 60-80%. Alternative centrality measures incor- 
porate forward neutral energy [|[ [l7|] and its correlation with 
multiplicity at midrapidity. The maximum variation of the p? 
spectrum for different centrality measures is 4% for central 
events and less than 4% for more peripheral events. This vari- 
ation is included in the systematic uncertainties of the reported 
spectra. 

Analysis of inclusive charged particle spectra for pt<2 
GeV/c has been described previously Accepted tracks 
for pt>2 GeV/c have |r;|<0.5 and a distance of closest ap- 
proach to the primary vertex less than 1 cm to reject back- 
ground. Acceptance, efficiency, and momentum resolution 
were determined by embedding simulated tracks into real 
raw data events. For pr>l-5 GeV/c, the Gaussian dis- 
tribution of track curvature k oc 1/pr has relative width 
dk/k = 0.016- (p T /(GeV/c)) + 0.012 for central events and 
5k /k = 0.011 • (pr/(GeV/c)) + 0.013 for peripheral events. 
Correction was made for distortion due to finite momentum 
resolution. 

The hadron yield decreases rapidly with increasing pt and 
its measurement is sensitive to small spatial distortions, which 
generate charge sign-dependent systematic errors in the mea- 
sured track curvature. Measurement of the summed hadron 
yield, (h + + h~)/2, is markedly less sensitive to such distor- 
tions than the yield of one charge sign alone. Each half of the 
cylindrical TPC is divided azimuthally into 12 sectors. High 
Pt tracks have small sagitta (s ~ 0.8 cm at pt=5 GeV/c) 
and are confined to a single sector. The sector-wise distribu- 
tion of (h + + h~)/2 has a pr -dependent rms variation of less 
than 5%, though with correlated variations for adjacent sec- 
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FIG. 1: Inclusive pr distributions of (h + + h~)/2. Non-central bins 
are scaled down by the indicated factors. The combined statistical 
and systematic errors are shown. Curves are fits to Eq. (j^). Hash 
marks at the top indicate bin boundaries for pr > 1 .5 GeV / c. 



tors. These effects contribute to the systematic uncertainty. 

The most significant background corrections are due to 
weak particle decays and anti-nucleon annihilation in detec- 
tor material. The former are estimated based on A, A and K® 
yields measured forpx < 2.5 GeV/c [ |T8[ [l9]] , with extrapo- 
lation to higher px using an exponential fit [ jl8| | . The latter are 
based on measured anti-proton yields [pop. 

The major correction factors and their uncertainties are 
given in Table |. "Tracking" incorporates efficiency, accep- 
tance, and the effects of the spatial non-uniformity of the TPC, 
with the latter dominating its systematic uncertainty. The to- 
tal systematic uncertainty of the spectra is the quadrature sum 
of the uncertainties in Table |. For the highest p? bins it is 
ss 27% for central events and sa 21% for peripheral events. 

Figure [l] shows the inclusive pr distributions of 
(h + + h~)/2 within |fy|<0.5 for various centrality bins, for 
Au+Au collisions at ^/snn=130 GeV. Error bars, which are 
dominated by systematic uncertainties at all pr, are generally 
smaller than the symbols. The spectra were fit by the pQCD- 
inspired power law function [ pT| 



dN 



2npT dpT 



c-(i+PT/p r\ 



(2) 



which describes px spectra of charged hadrons for NN col- 
lisions over a wide range of @ H ||]. Systematic 
changes in shape of the spectra with centrality are revealed 
by the fit parameters C, n, and mean transverse momentum 
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TABLE II: Geometric quantities, charged particle density per participant pair, and fit parameters to Eq. (g), for various centrality bins and for 
p±pat200 GeV[|l|, assuming <r™ e! =41 mb. The fits to the Au+Au data use uncorrelated measurement errors, which are largely systematic 
and non-Gaussian. Parameter errors shown also include correlated systematic uncertainties, which are added to parameter errors resulting from 
the fit. 













Power Law Fit (0.2 < p T < 6.0 GeV/c) 


centrality 


(N bin ) 


(Npart) 


dN ah 


2 dN ch 


C((GeV/c)- 2 ) 


<p T > (GeV/c) 


n 


dr) 


(Npart) d V 


0-5% 


965+ 67 


350I 4 


563±39 


3.22±0.23 


797±60 


0.520±0.010 


21.9±0.5 


5-10% 


764+59 


296+ 7 7 


452±32 


3.05±0.23 


654±50 


0.517±0.010 


20.7±0.4 


10-20% 


551+48 


232 +9 9 


344±24 


2.96±0.24 


520±40 


0.511±0.010 


18.9±0.4 


20-30% 




1651™ 


234±16 


2.83±0.26 


372±28 


0.504±0.011 


17.3±0.3 


30-40% 


210± 3 3 « 




144±10 


2.51±0.30 


252±19 


0.497±0.011 


17.2±0.3 


40-60% 


90±% 


621?! 


72.2±5.1 


2.35±0.42 


141±11 


0.480±0.011 


14.8±0.2 


60-80% 


20±l 


20±« 


21.0±1.5 


2.13±0.61 


50±4 


0.446±0.011 


13.0±0.2 


p + p(200) 


1 


2 


2.65 ±0.08 


2.65 ±0.08 


7.7±0.4 


0.392±0.003 


11.8±0.4 



<PT>=2po/(n — 3) in Table ||. Also shown in the Table 
are fit parameters for p ± p collisions at -^=200 Ge V J2l| ] . 
Eq. (Q) yields a poor fit relative to the uncertainties of the data 
for the most central bin, with systematic deviations from the fit 
function of 10-20%, whereas for more peripheral collisions it 
fits well, with parameters tending smoothly to those for p ± p 
collisions. 

A more direct comparison of yields for different centrali- 
ties relies on estimating (Nbin) and the mean number of par- 
ticipants (N part ) for each centrality bin. For this purpose, 
the distribution dcr/dNbin (and equivalently, d<r/dNp art ) was 
calculated using a Monte Carlo Glauber model Jl7[ ] with 
CT inei = 41 i 1 m b and Woods-Saxon nuclear matter density, 
using radius r = 6.5 ± 0.1 fm and surface diffuseness a = 
0.535 ± 0.027 fm. [Q. Percentile intervals of da/dN hin 
(calculated) and da/dN c h (measured) were equated to extract 
(Nbin) for each centrality bin. The Glauber model parame- 
ters and geometric cross section were varied to estimate the 
systematic uncertainties. Results are given in Table || 

Because charged particle distributions at midrapidity are 
used for centrality selection, biases in the relation between 
da/dN^ and da/dN c h due to fluctuations and autocorrela- 
tions were assessed. Variation of parameters in a Monte Carlo 
calculation of multiplicity fluctuations for fixed collision ge- 
ometry, and comparison of the measured multiplicity distribu- 
tion in an azimuthal quadrant centered on a high px particle 
against those in the other quadrants, both yielded negligible 
uncertainties in (Nbin) and (N part ). 

Table shows the charged particle yield per participant 
pair, ^ N 2 --T d ^ h , obtained by integrating fhepT spectra. The 
extrapolated yield in pr<0.2 GeV/c is ~ 20% of the total for 



all centralities. The dependence of 



dN, 



;n„ 



^ on (N part ) is 



consistent with observations in [|l7|]. 

Figure ^ shows the ratio of the central (0-5%) spectrum to 
that of the two peripheral bins (40-60%, 60-80%), normalized 
by (Nbin)- The dashed lines at unity and below show scaling 
with (Nbin) and (N par t) respectively, and the shaded regions 
show the systematic uncertainties from Table |l| The vertical 
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FIG. 2: Ratio of charged hadron yields within \r)\ < 0.5 for central 
over peripheral collisions, normalized to (Nbin). 



error bars on the data points are the uncertainties of the cen- 
tral data, while the horizontal caps are the quadrature sum of 
the uncertainties of both data sets. Approximate participant 
scaling at low px is seen. The ratio rises monotonically below 
Pt ~ 2 GeV/c and decreases at high px- The ratio of cen- 
tral over most peripheral achieves a value at pt=5.5 GeV/c of 
0.27±0.12 with additional uncertainty ±0.12 due to (N bi „), 
establishing significant suppression of charged hadron pro- 
duction at high pt in central collisions. 

Figure || shows Raa (pt) for various centrality bins relative 
to an NN reference spectrum, parameterized by Eq. (Q) with 

C°Zei = 2 ^-t mb/(GeV/c) 2 , Po = 1.9018$ GeV/c, and 



n = 12.98 



+0.92 
-0.47 



(the superscripts and subscripts are curves 



that bound the systematic uncertainty). The reference was de- 
termined by fitting Eq. (0) to UA1 p + p data at ^/s=200-900 
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FIG. 3: Raa{pt) for various centrality bins, for Au+Au relative 
to an NN reference spectrum. Error bars are described in the text. 
Errors between different pr and centrality bins are highly correlated. 



GeV [^TJ and extrapolating to ^=130 GeV. Extrapolation of 
the UA1 200 GeV spectrum to 130 GeV using pQCD calcula- 
tions agrees to within 5% with the reference at pt=6 GeV/c 
| p4| , ^5h , Correction to the NN reference for the UA1 accep- 
tance ( 1 77 1 < 2. 5), which differs from this analysis ( 1 77 1 < . 5), 
was based on two independent pQCD calculations p4| , |26| l, 
giving a multiplicative correction of 1.17 ± 0.06 at p^=2.0 
GeV/c and 1.35 ±0.10 at p T =5.5 GeV/c. The p T -dependent 
systematic uncertainty of the NN reference is the quadrature 
sum of the power law parameter and the acceptance correction 
uncertainties. Isospin effects are negligible forj»T<6 GeV/c 
[S], The error bars are the systematic uncertainties of the 
measured spectra, while the caps show their quadrature sum 
with the systematic uncertainty of the NN reference. 

Raa{pt) increases monotonically for j>t<2 GeV/c at all 
centralities, and saturates near unity for pr>2 GeV/c in the 
most peripheral bins. In contrast, Raa(pt) f° r me central 
bins reaches a maximum and then decreases strongly above 
Pt=2 GeV/c, showing suppression of the charged hadron 
yield relative to the NN reference of 0.36 ± 0.16 (sys) at 
Pt=5.5 GeV/c for the 0-5% bin, with additional uncertainty 
±0.03 due to (Nbin). Raa(pt) varies continuously as a func- 
tion of centrality, and no centrality threshold for the onset of 
suppression is observed [p7|]. 

In summary, charged hadron production in high energy col- 
lisions of heavy nuclei has been studied over a wide range 
of pt and event centrality. At high px, hadron yields scale 
with the number of binary collisions for peripheral collisions, 
while significant suppression of hadron production is seen for 
central collisions. This phenomenon indicates substantial en- 



ergy loss of the final state partons or their hadronic fragments 
in the medium generated by high energy nuclear collisions, 
though quantitative measurement of this effect requires addi- 
tional reference data. 
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